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Abstract 

Many extensions of the Standard Model predict super-weakly inter- 
acting particles, which typically have to decay before Big Bang Nu- 
cleosynthesis (BBN). The entropy produced in the decays may help 
to reconcile thermal leptogenesis and BBN in scenarios with gravitino 
dark matter, which is usually difficult due to late decays of the next- 
to-lightest supersymmetric particle (NLSP) spoiling the predictions of 
BBN. We study this possibility for a general neutralino NLSP. 

We elaborate general properties of the scenario and strong con- 
straints on the entropy-producing particle. As an example, we consider 
the saxion from the axion multiplet and show that, while enabling a so- 
lution of the strong CP problem, it can also produce a suitable amount 
of entropy. 



1 Introduction 



Thermal leptogenesis [1] is an attractive mechanism for generating the 
baryon asymmetry of the Universe, since it requires no additional ingre- 
dients beyond the see-saw scenario [2, 3, 4, 5, 6] introduced to explain the 
smallness of the observed neutrino masses. However, in supergravity theo- 
ries the required large reheating temperature results in a copious production 
of gravitinos [7, 8]. If the gravitino is heavier than other superparticles, this 
is problematic since it typically decays during or after Big Bang Nucleosyn- 
thesis (BBN) due to its extremely weak interactions. The energetic decay 
products can then cause unacceptable changes of the primordial light ele- 
ment abundances [9, 10]. If on the other hand the gravitino is the lightest 
supersymmetric particle (LSP), high temperatures tend to lead to a relic gra- 
vitino density exceeding the observed dark matter density [11]. In addition, 
the next-to-lightest supersymmetric particle (NLSP) becomes long-lived and 
thus can cause similar problems with BBN as an unstable gravitino [12]. 

Apart from abandoning supersymmetry or thermal leptogenesis, the gra- 
vitino problem has two further solutions. The first one is to reconcile thermal 
leptogenesis with a smaller reheating temperature. This is possible if there 
is a resonant enhancement of the generated asymmetry [13, 14, 15], some 
fine-tuning that violates the naturalness assumptions entering into the lower 
bound on the reheating temperature [16], or a violation of R-parity [17]. 

The second option is to accept a large reheating temperature and tackle 
the problems associated with the gravitino. One possibility is a very heavy 
gravitino that decays before BBN [18]. Alternatively, the gravitino could 
be very light, thus both avoiding a too large relic density and letting the 
NLSP decay before BBN [12]. The relic density of a gravitino LSP with 
mass in the GeV range can also be acceptable for rather large reheating 
temperatures [19], but in this case we have to protect BBN from the NLSP 
decays. This is possible if the NLSP decays relatively fast due to R-parity 
violation [20] or additional decay modes into hidden sector states [21, 22], if 
its decay products are only weakly interacting [23] or very low-energetic [24] 
and therefore harmless, or if its abundance is exceptionally small. The last 
alternative can occur for a stau NLSP [25, 26] in exceptional regions of 
parameter space, or for any NLSP whose abundance is diluted by entropy 
produced in late decays of another particle [27, 28, 29]. Of course, a combi- 
nation of the various options is possible [30, 31, 32, 33]. 

In this work, we study a scenario where the gravitino is the LSP and 
forms the dark matter. The abundance of the NLSP is decreased by late-time 
entropy production. This dilution has the additional motivation that long- 
lived particles are anyway present in many extensions of the Standard Model 
(SM) and usually have to decay before BBN to avoid problems, producing 
entropy in the process. The solution of the gravitino problem may thus 
be viewed as a complimentary by-product of such an extension rather than 
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an additional complication of the Minimal Supersymmetric Standard Model 
(MSSM). Of the several candidates for the NLSP, we consider the lightest 
neutralino. 

The paper is organised as follows. In Sec. 2, we review briefly the scenario 
of thermal leptogenesis with gravitino dark matter. Sec. 3 discusses the 
impact of entropy production on this scenario and the mechanism of entropy 
production by decaying matter. In Sec. 4, we present BBN constraints on 
a general neutralino NLSP after suitable dilution. We elaborate general 
properties of the scenario and strong constraints on the entropy-producing 
particle in Sec. 5, where we also consider as an example the saxion from the 
axion multiplet. 



2 Leptogenesis and Gravitino Dark Matter. . . 

In baryogenesis via standard thermal leptogenesis a cosmic lepton asym- 
metry is generated by CP-violating out-of-equilibrium decays of heavy 
right-handed Majorana neutrinos u^. Non-perturbative sphaleron pro- 
cesses [34, 35] convert the lepton asymmetry into a baryon asymmetry t/b- 
In the case of hierarchical masses the maximal resulting baryon-to-photon 
ratio of the Universe can be given as [36, 37] 

for the MSSM and weak washout, while the observed value lies in the range 
5.89 x 1(T 10 < r? B bs < 6.49 x 1(T 10 (2a) [38]. M v i denotes the Majorana 
mass of the lightest right-handed neutrino. The given baryon asymmetry is 
maximal in the sense that the CP violation in the decays is chosen to be 
maximal [39]. Since thermal leptogenesis strongly favours hierarchical light 
neutrino masses [40] , the mass m„3 of the heaviest left-handed neutrino has 

to be close to \J Am^ ~ 0.050 eV, using the best-fit value from neutrino 
data [41] and assuming a normal mass ordering. The efficiency factor kq 
should be computed case-by-case by solving the relevant Boltzmann equa- 
tions [42, 43, 44]. For zero initial abundance in the small M u i regime [36], 

R. 

i.e. for M u i < 4 x 10 13 GeV, the maximal value is «g eak ~ 0.18 [45]. This 
value is reached for 

_ 87r 2 ,/ft,r v 2 n 

m , A ~ = P ~ 1.6 x 10" 3 eV , 2) 

3-/L0 Mpi K J 



where m* is known as the equilibrium neutrino mass, v ~ 174 GeV, and 
Mpi ~ 2.44 x 10 18 GeV. The effective neutrino mass 

( m D m D)n 

m ,A = (3) 

l M..\ y 1 

R 
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equals the mass of the lightest neutrino if the Dirac mass matrix m-rj is 
diagonal. Its natural range is m„i < m„i < m u s. The parameter A denotes 
the dilution factor by entropy production after the decay of the right-handed 
neutrinos. It equals one in standard cosmology, while we will consider the 
general case A > 1 later on. 

There are some uncertainties entering (1). Possible spectator field uncer- 
tainties [46] and flavour effects [47, 48] are neglected, and the naive sphaleron 
conversion factor [49, 50] is used. We have assumed the particle content of 
the MSSM with = 228.75 for the number of effectively massless degrees 
of freedom at high temperatures. To be conservative we consider the effects 
of the MSSM by a factor 2\/2 relative to the SM, which is valid for weak 
washout [37] . For strong washout this factor reduces to y/2. 

We see from (1) that leptogenesis in its minimal version as described 
above can generate the observed baryon-to-photon ratio of the Universe, 
because r/g 1 " can exceed r/g bs . On the other hand, it is clear that there is a 
lower bound M v x > 2 x 10 9 GeV. 

It is especially appealing for the considered neutrino mass range that lep- 
togenesis can emerge as the unique source of the cosmological baryon asym- 
metry [40]. Wash-out processes may reduce a pre-existing asymmetry by 
two to three orders of magnitude for the situation of (1). Stronger washout 
decreases the efficiency factor and thus requires a larger right-handed neu- 
trino mass to keep 77g ax > r/g bs . For thermal leptogenesis, the bound on the 
lightest right-handed neutrino mass can be translated into a lower bound on 
the reheating temperature after inflation, Tr > M u i . In the strong washout 
regime, i.e. for fh u i > m*, this changes to Tr > 0.1 M v i [45], but we cannot 
relax the bound on the absolute value of Tr, since in this case the efficiency 
factor decreases as well, requiring a larger M v \ . 

The required high temperatures also lead to thermal production of a 
significant gravitino relic density 

nf /2 h* = m 3/2 Y$ 2 (T ) S P^, (4) 

where s(To) refers to today's entropy density of the Universe. Together 
with the Hubble constant h in units of 100 km Mpc -1 s -1 and today's crit- 
ical density po, we obtain s(To)h 2 / po ~ 2.8 x 10 s GeV -1 . The gravitino 
abundance for low temperatures T\ ow <C 7r is given by [51, 52] 




where the gauge couplings g\ = {g',g,g s ), the gaugino mass parameters 
Mi as well as the constants k { = (1.266,1.312,1.271) and yi/10- 12 = 
(0.653,1.604,4.276) are associated with the gauge groups U(1)y, SU(2)l 
and SU(3)c, respectively. 
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Without entropy production, the gravitino yield from thermal produc- 
tion at the present temperature Y^? 2 (To) = Fg? 2 (Ti ow ). With entropy pro- 
duction after the gravitino production in the early Universe, 

y 3 t / p 2 (To) = A- 1 y 3 t ; 2 (r low ). (6) 

As mentioned before, A = 1 in standard cosmology. For reasons that we 
will explain below, we will consider late-time entropy production at T <C 
M u i , T\ ow . Then the same dilution factor A appears in (1) and (6). 1 

R. 

From (4) and (5) we see that since the gravitino is the LSP, for fixed 
gaugino masses the relic gravitino density typically decreases for increasing 
gravitino mass. Assuming universal gaugino masses at the GUT scale, we 
can approximate 

0* ~ 11 A- 1 ( *a ( M 9^) \ 2 (WGeV) 

^ - U - iiA V3x 108 GeVy* VlO 3 GeVy* V ™3/2 )' { > 

Thus, for given reheating temperature and gaugino masses, we obtain a lower 
bound on the gravitino mass exploiting the requirement tl^^h 2 < firjM^ 2 = 
0.112 ±0.007 (2a) [38]. 

In order to summarise the issues discussed so far, we combine (1) and 
(7) using the best-case relation Tr ~ M u i to eliminate the right-handed 
neutrino mass, arriving at 




Note that the dilution factor A cancels out. Recalling the discussion af- 
ter (1), m u 3 cannot be raised without lowering kq. Thus, even for the most 

optimistic scenario with Tr = 2 x 10 9 GeV the gravitino mass is restricted to 
a rather large value > 40 GeV. In other words, there is considerable tension 
between thermal leptogenesis and gravitino dark matter. 

Even worse, the NLSP decay problem is a definite clash between both 
notions, since gravitino LSP masses larger than about 10 GeV are excluded 
in most cases. In the MSSM with conserved R-parity, the NLSP has to decay 
into the gravitino and SM particles. It decays typically with a long lifetime 
due to the extremely weak interactions of the gravitino. If these decays occur 
during or after BBN, the emitted SM particles can change the primordial 
abundances of the light elements [12, 53, 54]. Specific setups like the Con- 
strained MSSM have been studied in [55, 56, 57, 58, 59]. Only in exceptional 

1 For Tr M„i , it could be possible to produce entropy in between, so that only 
would be diluted But not t)b- 
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regions of the parameter space [24, 25, 26], a stau NLSP not orders of mag- 
nitude heavier than the gravitino could be consistent with all constraints, 
allowing reheating temperatures Tr ~ 10 9 GeV. Generically a conservative 
upper bound on the reheating temperature Tr < few 10 8 GeV is found. 
Since the maximally produced baryon asymmetry is too small, various NLSP 
candidates have been investigated more model-independently [60, 54] to 
identify best-case scenarios: a sneutrino [23, 61, 62, 63, 64, 65] actually 
could allow large enough reheating temperatures with reasonable masses 
due to its invisible decays. A stop [66, 67] or a general neutralino [68] are 
not reconcilable with thermal leptogenesis for masses below a TeV. They 
would definitely require m 3/ / 2 < 10 GeV. 

Altogether, the strongest conflict between thermal leptogenesis and gra- 
vitino dark matter is found in the NLSP decay problem. This is embodied 
in (8) by the restriction to small gravitino masses, m 3 / 2 < 10 GeV. 

Entropy production after the freeze-out of the NLSP dilutes its density. 
Thus, late-time entropy production can naively resolve this conflict for any 
NLSP within or without a specific model. The relic density prior to its 
decay, 

^ nlsp = A- 1 ^° lsp ' ( 9 ) 

is reduced compared to its freeze-out density f^ sp by ^ ne dilution factor A, 
which is the same as in (1) and (6). In Sec. 4 we show how BBN constraints 
on a general neutralino with a gravitino LSP with m 3 / 2 = 100 GeV are 
softened by A > 1. 

3 ... with late-time Entropy Production 

3.1 Thermal Leptogenesis and Gravitino Yield with late- 
time Entropy Production 

From (1) we see that a significant dilution A > 1 can only be compensated 
by a larger M v i , since all the other parameters are chosen already to be 

R, 

optimal. Due to the requirement Tr > M v i this gives a linear shift of the 
required reheating temperature. Since the gravitino density depends also 
linearly on the reheating temperature (5) and is diluted in the same way as 
the baryon density, such a compensation seems to give a trivial shift of the 
problem to higher reheating temperatures. However, there are aspects that 
do not show up in (1) and (8), in addition to the impact on the NLSP decay 
problem. 

Most importantly, in the domain of large M u i washout processes reduce 
the efficiency factor kq exponentially. In the case of hierarchical neutrinos, 
this domain corresponds to M v i > 4 x 10 13 GeV. From this we would obtain 

R 

A < 2 x 10 4 . However, while at low M v \^ many numeric examples and an 
analytic approximation for kq [45] exist in the literature, the situation for 
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larger M v i is less well-studied. As an additional complication, for T S> 

10 9 GeV more spectator processes are in equilibrium and thus should be 
taken into account. Hence, it is not clear whether the maximal value KQ 6ak 
can be reached for M u i ~ 4 x 10 13 GeV. Consequently, we expect an upper 
bound 

A < A max ~ 10 3 . . . 10 4 . (10) 

We would like to stress that this is an intrinsic bound of the problem. It is 
stronger than bounds from perturbativity of Yukawa couplings (A < 10 5 ) or 
the requirement of a reheating temperature below the GUT scale (A < 10 7 ). 

We remark that according to Fig. 6b of [36] there is a much stronger 
bound with roughly A max < 10 2 for quasi-degenerate neutrino masses. 
Thus, thermal leptogenesis with late-time entropy production requires hier- 
archical neutrinos even more than thermal leptogenesis already does. 

Late-time entropy production leads to a strong reduction of the allowed 
parameter space for successful thermal leptogenesis. Since the required min- 
imal M v i is increased, the range of allowed values for kq and the neutrino 
mass parameters is reduced. However, the same region of parameter space 
is already favoured by the need to keep the reheating temperature as low as 
possible in order to avoid the overproduction of gravitinos. Therefore, late- 
time entropy production does not reduce the parameter space of thermal 
leptogenesis with gravitino dark matter. 

In (5) one has to consider the impact of the running couplings and masses 
due to the shift of the reheating temperature. For example, if we increase 
Tr from 3 x 10 9 GeV to 3 x 10 13 GeV and choose A = 10 4 to compensate, 
f2 3 ^ 2 decreases by 25%. 2 Note that this effect is unavoidable and softens 
the tension between thermal leptogenesis and gravitino dark matter already 
before considering the impact of entropy production on the NLSP decay 
problem. 

Another possibility is a gravitino with such a small mass that it comes 
into thermal equilibrium after reheating. Then its relic abundance becomes 
independent of the reheating temperature, which allows Tr 3> M u i . Taking 
into account the lower limit on the mass of a warm dark matter particle [69] , 
it turns out that its relic energy density exceeds the observed dark matter 
density in standard cosmology. However, already a A of a few dilutes the 
gravitino sufficiently to make it viable warm dark matter again [31, 32]. For 
A ~ 10 3 it forms cold dark matter with m 3 /2 — 1 MeV [33, 70, 71]. Note 
that for these small masses the NLSP decays before BBN, so that the decay 
problem is absent. 

2 Besides, the electroweak contributions double their contribution to the total yield to 
about 30%. 
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3.2 Entropy Production by decaying Matter 



In this section we discuss briefly how decaying matter can produce consid- 
erable entropy in the early Universe [72, 73]. We consider a non-relativistic 
and long-lived particle species <j) with chemical potential p = in a flat 
Friedmann- Robertson- Walker Universe. When <fi drops out of chemical equi- 
librium, its abundance = n^j ' s "freezes out" , where n$ denotes its number 
density and 



flaton decay or thermally after reheating, if (f> never enters chemical equi- 
librium. The contribution of non-relativistic particles to the energy density 
decreases as p ma t oc R~ 3 , where R denotes the scale factor. Since the energy 
density of radiation in the Universe, 



decreases oc i? -4 , Pmat/Vrad grows oc R. Since R grows with time, at some 
time or temperature the unstable species (j) comes to dominate the 
energy density automatically, if its lifetime > £T. If the Universe has 
been dominated by radiation before, it enters a phase of matter domina- 
tion that lasts roughly till (p exponentially decays at r^. Here we assume 
that everything released by the particle decay is rapidly thermalized, i.e. 
on timescales At <C H^ 1 ~ expansion time. At some intermediate time 
* — t$ i T <t>/^ ) 3 ^ 5 ^ ne radiation produced in decays of cf> starts to become the 
dominant component of the radiation energy density. The temperature of 
the Universe begins to fall more slowly, T oc i?~ 3 / 8 , than the usual T oc Rr 1 . 
The Universe is never reheated, since the temperature decreases at all times. 
From t ~ (f^/t^) 3 ^ till t ~ r^, the entropy per comoving volume 5 is 
growing oc i? 15 / 8 . At the Universe becomes purely radiation-dominated 
again with T oc R^ 1 and a temperature T^ ec = T(T^,)\ Tad _ dom , wnere we use 
the time-temperature relation for a radiation-dominated Universe, 



This is the temperature after significant entropy production T atter ; which 
would be identified as the reheating temperature in the approximation of 
simultaneous decay of all cf> particles. We identify T after = T^ ec . If never 
dominates over p ra d, </> decays never produce a significant amount of entropy 
relative to the initial entropy. Then the produced entropy is negligible. 

3 For simplicity we use g* only, since the temperatures occurring in this work are above 
1 MeV, where g*s = g*- 




Prad = — s*(T)T 4 



(12) 




(13) 
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However, if p^ dominates over p ra 4 before r^, the produced entropy dilutes 
significantly any relic density by a factor A. 

The dilution factor A is defined as the ratio of entropy per comoving 
volume after <j) decay Sf over the initial entropy per comoving volume Si 
and can be expressed as 

A = ^.0.82< 5 yy/ 4 ^f!, (14) 

1/3 

where the angle brackets indicate the appropriately-averaged value of gj 
over the decay interval. We see how A is determined by the properties of the 
unstable particle, i.e. its mass and lifetime ta. Meanwhile it is assumed 
that 70 > t^. The pre-decay abundance Ya of the unstable particle depends 
on both its interactions and the earlier cosmology. 

For convenience we would like to rephrase (14) in terms of tempera- 
tures. Without entropy production after the generation of the pre-decay 
abundance, it is constant till the particle decays. With p = ran we find 



Prad 



<t> 



where we have used (11), (12) and p<j, = p va( \ at . 

Using (13) we can replace the particle lifetime in (14) as 



45 

27^ 



V gJiT^M^T^y 1 . (16) 



Plugging (15) and (16) into (14) we obtain 

/ 1/3 \ 3/4 rp= 

A = 0.75 -f^ — (17) 

y* K 1 ^ i 4> 

This linear growth in temperature can also be expressed in terms of energy 
densities, since 



n^rriA = sY^m^ = ^— - g*(T) T 3 T, , (18) 



2vr2 3 
45 4 

where we have used (11) and (15). Taking this together with (12) we find 

Da, T, 



Pra,d T 

Thus, for T = Tf ec we see that 



(19) 



^ c =jHn ec ), (20) 
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where p ra( j is the density of the "old" radiation, i.e. it does not include the 
radiation from cp decays. 

The standard Big Bang model has been tested thoroughly up to tem- 
peratures around 1 MeV, where BBN occurs. Investigations of the ther- 
malization of neutrinos produced in (ft decays or subsequent thermalization 
processes lead to lower limits on the temperature of the Universe after the 
entropy production T a [74]. Neutrinos, which can thermalize through 
weak interactions only, are most important. All other SM particles thermal- 
ize much faster due to their stronger interactions. The bounds found are in 
the range 

T$ ec > T*f n er ~ (0.7 ... 4) MeV (21) 

where weaker bounds come from BBN calculations [75, 76] and stronger 
bounds rely on the neutrino energy density [77, 78] exploiting overall best 
fits for cosmological parameters. We take T^ ec > 4 MeV ~ Tbbn as lower 
bound. 

Going back in time, thus towards higher temperatures, the first cosmo- 
logical event important for our scenario is the freeze-out of the NLSP. 4 Stan- 
dard computations of relic abundances rely on the assumption of radiation 
domination during freeze-out. If the Universe is dominated by matter during 
NLSP freeze-out, the NLSP relic abundance is increased. Taking the later 
dilution by entropy production into account, the overall effect remains a re- 
duction [79] . The effects of different cosmological scenarios on relic densities 
have been studied [80, 81, 82] and there are computer codes [83]. In partic- 
ular, the neutralino has been investigated, also considering the production 
of neutralinos in the decay of a dominating matter particle [84, 85, 86, 87]. 
Since it is the easiest case to study, we take < T^ sp . Thereby the 
Universe is radiation-dominated during NLSP freeze-out happening at T^° sp 
and the standard computations hold. 5 Later we will find that the window 
between BBN and NLSP freeze-out is favoured intrinsically by the scenario. 

Sticking to this particular window, we can evaluate (17), 

A * 0.75 x 10 3 ( ( ^) > (22) 

VlOO GeV/ I T^ cc J y ' 

where we have plugged in = T^ sp ~ m n \ sp /25 and (g¥ 3 ) — 2.2 ~ 

gl /3 {T$ cc ) with g*{T$ ec ) = 10.75, exploiting the fact that for 4 MeV <T< 
4 GeV the effective relativistic degrees of freedom are known [88]. If we 

4 The QCD phase transition occuring between BBN and NLSP freeze-out seems not 
to deliver any constraint on our scenario. 

5 Using the simple estimate H(T^° sp ) ~ r(Tn° sp ), where T is the rate of NLSP anni- 
hilations, one finds that for = T^ Bp the NLSP abundance is increased by a factor 
of only v2 compared to the standard case of radiation domination, while the freeze-out 
temperature stays nearly constant. 
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compare (22) and (10), we see that the cosmological window between BBN 
and NLSP freeze-out is not only the first and easiest but also sufficiently 
large to produce enough entropy to come close to the upper limit on A set 
by thermal leptogenesis itself. 

This discussion assumes that there is no further entropy production after 
the generation of Y^. Otherwise, would be diluted like any other relic 
abundance, i.e. Y^ — >■ Y[ = A± 1 Y ( j > . There are two possibilities for the 
impact of such an earlier entropy increase Ai > 1. i) Despite the dilution, 
<p dominates the Universe for some time. Then the later entropy production 
by the decay of 4> is simply reduced by a factor Ai, as we see from (14) 
since lifetime and mass of the unstable particle are unchanged, ii) The 
relic abundance of cfi becomes so small that the particle never dominates the 
energy density of the Universe. Then (14) does not hold and S — const., i.e. 
A = 1. 

After an arbitrary number of late events of entropy production Aj la- 
beled by i = 1, 2, . . ., where the index implies a time-ordering with larger i 
corresponding to later decays, the total dilution factor is 

Atot = n A »- ( 23 ) 

i 

Here, "late" indicates that all decays happen after the freeze-out of all un- 
stable particles supposed to produce significant entropy, so that their relic 
abundances are diluted by each earlier decay. This implies 

where Aj(A,-<j = 1) refers to the dilution factor obtained from (14) without 
considering the other dilutions in the calculation of Y^. As mentioned, we set 
Aj = 1, if a decaying particle does not come to dominate the energy density 
of the Universe. One can easily convince oneself that the total dilution is 
simply given by the largest individual dilution factor, 

A tot = maxjA^A^i = 1)} . (25) 

The upper bound (10) limits Atot- The dilution of the NLSP abundance 
can be smaller than A to t if some decays happen before NLSP freeze-out. 
Thus, we see from (24) with (10) how our requirement of sufficient entropy 
production after NLSP freeze-out restricts the possibility of earlier entropy 
production. 

4 BBN Constraints on a Diluted Neutralino NLSP 

In this section we present constraints from Big Bang Nucleosynthesis on a 
neutralino NLSP in the case of a gravitino with a mass of m.3/2 = 100 GeV 
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being the LSP. We investigated those bounds excluding the possibility of 
entropy production in [68] and found for masses below a TeV a maximal 
gravitino mass of a few GeV. In the following we assume that the neutralino 
is diluted after its freeze-out by a factor A = 10 3 . It is trivial to infer the 
impact of arbitrary As. BBN constraints on a stau NLSP with A up to 
2 x 10 have been studied in [27, 28], where it has been found that interesting 
parameter regions are allowed for dilution factors A ~ 10 3 . 

For a discussion of the neutralino decay channels, branching ratios and 
more details we refer to [68]. To determine model-independent constraints 
within the MSSM we take all points that are not ruled out by LEP up to a 
mass of 2 TeV, while we fix the masses of the sfermions to be above 2 TeV. 
To keep our analysis as general as possible we do not fix all supersymmet- 
ric parameters according to a specific scenario, but instead we set the soft 
SUSY breaking parameters at the low energy scale. We keep the majority of 
the parameters fixed and vary the gaugino and Higgsino mass parameters to 
study how the lifetime and number density vary with the mass and composi- 
tion of the lightest neutralino. We plot these points against the hadronic and 
electromagnetic BBN bounds in Figs. 1-3. The bounds are taken from [53] 
and the different curves are explained in the figure caption. The vertical 
axis corresponds to the fraction of the energy density that decays to electro- 
magnetic or hadronic products. A A > 1 shifts all points downwards on this 
axis by a factor of A. Therefore it is easy to infer constraints for arbitrary 
A once a plot with fixed A is given. 

Hadronic bounds are generally more constraining. However, it has been 
found that large gravitino masses, for which light neutralinos have a low 
hadronic branching ratio, are excluded by the electromagnetic bounds. 

In Fig. 1 we consider a mixed bino-wino NLSP. The large dip corre- 
sponds to resonant annihilation into the pseudo-scalar Higgs, which hap- 
pens for our choice of parameters at a neutralino mass m x ~ 1150 GeV. To 
increase r^g ax we are more interested in the region of small NLSP masses, 
since small NLSP masses allow more easily for small gluino masses in (8). 

Thanks to the dilution by entropy production the wino overcomes the 
electromagnetic bounds for any mass even for masses close to the gravitino 
mass. If the neutralino is mainly bino the electromagnetic bounds are more 
involved. For a bino-like neutralino, masses below about 450 GeV are ex- 
cluded. Smaller and smaller masses become allowed when the wino compo- 
nent increases, so that there is allowed space for binos with a non-negligible 
wino component and m x ~ 200 GeV or even smaller masses. 

The hadronic bounds exclude most of the parameter space for a bino- 
wino with dominant bino component even with A = 10 3 . The mixed bino- 
wino states with m x ~ 200 GeV mentioned above are found on the less 
conservative 6 Li/ 7 Li exclusion line for a decaying particle of 100 GeV mass. 
Thus we find many points that should not be considered as strictly ex- 
cluded with masses around 200 GeV and also mixed bino-wino states that 
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Hadronic BBN bounds 
m3, 2 = 10OGeV with: M 3 = 2200, \i = 2200, taig = 10, agn(|j)= 1. 




-2 -1 1 2 3 4 5 6 7 8 9 10 11 12 -2 - 1 1 2 3 4 5 6 7 8 9 10 11 12 

lcg 10 (T,,/sa:) log, (T x /sec) 



Figure 1: Energy density of the bino-wino neutralino decaying into electro- 
magnetic/hadronic products compared with the BBN electromagnetic (left) and 
hadronic constraints (right) for the case of a 100 GeV gravitino mass and a dilu- 
tion factor A = 10 3 . The bounds are taken from [53]: the continuous (dashed) 
lines correspond to more (less) conservative bounds for the 6 Li to 7 Li ratio, and 
the region between the curves should not be considered as strictly excluded. The 
red/upper and violet/lower curves in the hadronic plots are the constraints for 
1 TeV and 100 GeV decaying particle mass, respectively. The mass increases from 
right to left as heavier particles decay faster. The composition goes from bino at 
the top to wino at the bottom while the colours give the dominant component. 
The deformation between the left and right panel is due to the mass dependence of 
the hadronic branching ratio with lighter NLSPs having lower branching ratios to 
hadrons. In contrast the electromagnetic branching ratio is always nearly one. 

are allowed with masses smaller than 200 GeV. Winos with m v < 400 GeV 
overcome even any less conservative bound. For 400 GeV < m x < 1100 GeV 
the wino could violate the less conservative bound, while even larger masses 
become allowed again. 

Altogether, the situation is qualitatively different for bino and wino. 
While the wino safely overcomes all bounds, especially at low masses, a bino- 
like neutralino with reasonable mass stays excluded even for much larger 
dilution factors that would be in contradiction with successful thermal lep- 
togenesis (10). However, there is also some space for bino-wino mixed states 
that are mainly bino with masses below 200 GeV. 

In Fig. 2 we consider a mixed bino-Higgsino NLSP. The dip is broader 
in this case, since the Higgsino component that couples to the pseudo-scalar 
Higgs is larger. Thanks to the dilution the Higgsino overcomes the electro- 
magnetic bounds like the wino for all masses, even though A should not be 
much smaller than roughly 10 2 to allow for light Higgsino neutralinos. For 
the bino the situation is comparable to the case of mixed bino-wino. No 
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Figure 2: Energy density of the bino-Higgsino neutralino decaying into electro- 
magnetic/hadronic products compared with the BBN electromagnetic (left) and 
hadronic constraints (right) for the case of a 100 GeV gravitino mass and a dilu- 
tion factor A = 10 3 . The bounds are taken from [53]: the continuous (dashed) 
lines correspond to more (less) conservative bounds for the 6 Li to 7 Li ratio, and 
the region between the curves should not be considered as strictly excluded. The 
red/upper and violet/lower curves in the hadronic plots are the constraints for 
1 TeV and 100 GeV decaying particle mass, respectively. The mass increases from 
right to left as heavier particles decay faster. The composition goes from bino at 
the top to Higgsino at the bottom while the colours give the dominant component. 
The deformation between the left and right panel is due to the mass dependence of 
the hadronic branching ratio with lighter NLSPs having lower branching ratios to 
hadrons. In contrast the electromagnetic branching ratio is always nearly one. 

mixed bino-Higgsino state with a dominant bino component is allowed with 
masses as low as 200 GeV, though. 

Again, the hadronic bounds exclude most of the bino parameter space. 
Exceptions are found in the dip and at very large masses. There are states 
with comparable bino and Higgsino components and m x > 200 GeV — thus 
not excluded by the electromagnetic bounds — violating the less conservative 
hadronic bound. Higgsino neutralinos lighter than 250 GeV escape even 
these constraints, while they are excluded for 670 GeV < m x < 1100 GeV. 

Altogether, we find that for mixed bino-Higgsino only states that are 
mainly Higgsino allow for preferrable small masses but then even down to 
the gravitino mass. Considering only the conservative hadronic bound from 
the 6 Li to 7 Li ratio, in addition maximally mixed states with masses in the 
region around 230 GeV become allowed. A > 10 3 would allow for larger 
bino components in the mixed bino-Higgsino. 

In Fig. 3 we consider a mixed wino-Higgsino NLSP. Thanks to the 
dilution it overcomes the electromagnetic bounds for all masses and mixings. 
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Figure 3: Energy density of the wino-Higgsino neutralino decaying into electro- 
magnetic/hadronic products compared with the BBN electromagnetic (left) and 
hadronic constraints (right) for the case of a 100 GeV gravitino mass and a dilu- 
tion factor A = 10 3 . The bounds are taken from [53]: the continuous (dashed) 
lines correspond to more (less) conservative bounds for the 6 Li to 7 Li ratio, and 
the region between the curves should not be considered as strictly excluded. The 
red/upper and violet/lower curves in the hadronic plots are the constraints for 
1 TeV and 100 GeV decaying particle mass, respectively. The mass increases from 
right to left as heavier particles decay faster. The composition goes from Higgsino 
at the top to wino at the bottom while the colours give the dominant component. 
The deformation between the left and right panel is due to the mass dependence of 
the hadronic branching ratio with lighter NLSPs having lower branching ratios to 
hadrons. In contrast the electromagnetic branching ratio is always nearly one. 



We are especially interested in the small mass region. By vertical shifts of 
all points we can search for the minimal dilution factor A mm to overcome 
the electromagnetic bounds at small masses. A wino close to the gravitino 
mass becomes allowed for 

Af * ~ 25 . (26) 
Larger A allows for more wino masses and eventually for light Higgsinos at 

Af n ~ 90 . (27) 

Considering the hadronic constraints, winos with m x < 400 GeV and 
Higgsinos with m x < 200 GeV satisfy all bounds. Wino-Higgsinos with 
larger mass can be in conflict with the less conservative bound, and the 
wino overcomes the more conservative one completely. Disregarding the 
dip, Higgsinos are excluded for a window 700 GeV < m x < 1300 GeV. 

In summary, entropy production after NLSP freeze-out can allow for a 
gravitino LSP of 100 GeV mass with a light neutralino NLSP. This recon- 
ciles thermal leptogenesis and gravitino dark matter within the scenario of 
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a light neutralino NLSP. However, this depends on the composition of the 
lightest neutralino. The wino is the best case due to its small freeze-out 
abundance. Also a light Higgsino becomes allowed for reasonable dilution 
factors. A light bino-like neutralino, which is typical for the Constrained 
MSSM, stays excluded for m 3 / 2 = 100 GeV even if the possibility of entropy 
production after freeze-out is exploited. 

Coming back to thermal leptogenesis, strictly speaking none of the points 
in the figures allows for a sufficiently high reheating temperature, since the 
gluino mass has been fixed at 2.2 TeV. However, in the parameter space re- 
gions with smaller neutralino masses, the gluino mass can be lowered without 
affecting our considerations at all. Consequently, all allowed points with a 
neutralino mass below a TeV can be compatible with thermal leptogenesis. 

5 Search for a viable Candidate 
5.1 General Requirements 

In this section we discuss candidates for the entropy-producing particle 4> of 
the previous sections. After enumerating the required properties in general, 
we exemplify in detail an implementation of the scenario with the axion 
multiplet. 

To dilute the NLSP relic density, <p must i) decay after NLSP freeze-out. 
But, for sure, it ii) decays before BBN. Thus the lifetime or equivalently 
the decay temperature T^ ec is constrained to a window. The particle has 
to be Hi) produced in the early Universe such that it dominates the energy 
density before BBN. Meanwhile we stick to the case where its relic density 

iv) does not come to dominate before NLSP freeze-out. Thus the relic density 
prior to its decay = Y^m^s is also constrained to a window. Requirements 
iii)+iv) imply that the dominance of 4> has to grow with the expansion, which 
is true for non-relativistic matter. So 4> is implicitly assumed to become non- 
relativistic before BBN. The requirements i)+h) and iii)+iv) constrain two 
different quantities and p<p, which are determined in different ways but 
by the same properties of (p, namely its couplings and mass. 

Requirements i)+h) constrain only the total decay rate r^ ot = tT 1 . In 
fact, the branching ratios of (ft into the LSP and NLSP are also constrained. 
The branching ratio into the NLSP i?0^ n i sp+ ... must be so small that the 

v) NLSP decay problem is not reintroduced by the decay. Branching ratios 
into the LSP are always restricted by overproduction. Especially when the 
LSP is already produced thermally as in our scenario of thermal leptogenesis 
with gravitino dark matter, (ft should vi) not produce too many LSPs in its 
decays. Since <f> even dominates the energy density of the Universe at its 
decay, the requirements v)+vi) force the corresponding branching ratios to 
be — at least — close to zero. 

In addition, <p must be vii) compatible with gravitino dark matter. For 
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No. 


Requirement 


Comment 


i 


Txlsc ^ r f~ , IO 
1 4> < J nlsp 


to nave eitect on ii n i sp 


ii 


rnn f^c m 

J-tj, > -tBBN 


1 J ■ ~t I A T "A TV T 

not to spoil BBN 


iii 


P<t> (rpAGC\ v. 1 
Prad^ ^ > 1 


O(10) < A < 10 4 


iv 


Prad UnlspJ < 1 


for standard NLSP freeze-out 


V 


£^->nlsp+... — 


from NLSP decay problem 


vi 




from overproduction (0 3 y 2 ~ SIdm) 


vii 


e-g- r 3/2 > t 


compatibility with gravitino dark matter 


viii 


ii) and v)— vii) 


for by-products; no new problems 



Table 1: List of requirements for our scenario of entropy produced by <f> to dilute 
the NLSP. 

example, the gravitino would become unstable due to the existence of (ft, if it 
could decay into (ft. This would take away the explanation for the observed 
dark matter abundance, if the gravitino lifetime were too short, or by itself 
be in conflict with other observations. 

Finally, viii) unavoidable by-products of (ft have to be harmless. For ex- 
ample, such by-products are the supermultiplet partners in SUSY. They are 
harmless, if they do not violate ii) or vii), are free of the problems solved by 
v)+vi) and do not introduce new problems on their own. 

Altogether, the properties of (ft seem to be highly constrained. We sum- 
marise the requirements in Tab. 1. The number of free parameters — mass 
and couplings — is finite. Since they enter in different ways for different con- 
strained quantities, it is not a matter of course that the scenario of late-time 
entropy production is viable at all. Especially if is produced thermally 
via scatterings, the same coupling might be responsible for the production 
and the late decay. 

On the other hand, many extensions of the SM contain or predict super- 
weakly interacting and hence long-lived particles. Such particles generically 
satisfy i), if not by definition. In order to ensure that they are harmless, 
one usually demands that they decay before BBN conform to ii). Thermal 
leptogenesis places the upper limit (10) on the maximally allowed dilution, 
which implies that for decay right before BBN iv) has to hold at least ap- 
proximately. Considering high reheating temperatures and the growth of 
Pmat/Vrad °^ R> if is probable that the energy density of late-decaying par- 
ticles dominates over the radiation energy density at their decay. Thus, 
iii) can be considered as fulfilled generically, which in fact normally poses a 
problem. Besides, the decay into superparticles usually has to be suppressed 
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in order to avoid producing too much dark matter and further late-decaying 
particles like the NLSP or the gravitino, in case it is not the LSP. Conse- 
quently, v) and vi) are generic, too, possibly amended by -B«^i sp +.. . - 0, if 
the gravitino is not the LSP. In any case the scenario has to be compatible 
with whatever is supposed to form the dark matter, so that vii) is generic. 
Also viii) arises as a generic requirement on any late-decaying particle and 
is particularly constraining in supersymmetric models. 

In summary, <fi is severely constrained such that the scenario of entropy 
production to dilute the NLSP density might appear unappealing. How- 
ever, in extensions of the SM containing long-lived particles, in principle 
i)+ii) and v)-viii) are no new requirements and are present — in appropri- 
ate form — without considering entropy production at all. If T^ cc ~ Tbbn! 
successful thermal leptogenesis favours the situation of iv). Finally, for the 
corresponding high reheating temperatures hi) is generic. Thus, all the re- 
quirements of Tab. 1 either have to be fulfilled or are generically fulfilled. In 
other words, the solution of the generic problems of long-lived particles may 
well cause the entropy production desired to solve the NLSP decay problem 
and thereby reconcile thermal leptogenesis and gravitino dark matter. 

With a specific candidate at hand the details have to be worked out. One 
has to determine whether a candidate is excluded, not useful or can be the 
solution and how generically this is true. As an example, we investigate the 
axion multiplet, which is motivated by a completely disconnected problem 
of the SM. 

5.2 Example: Axion Multiplet 

The strong CP problem of the SM can be solved by the Peccei-Quinn (PQ) 
mechanism [89, 90]. An additional global U(1)pq symmetry referred to as 
PQ symmetry, which is broken spontaneously at some PQ scale, can explain 
the smallness of the CP-violating O-term in QCD. The pseudo Nambu- 
Goldstone boson associated with this spontaneous symmetry breaking is 
called axion [91, 92]. It has not yet been observed. However, axion physics 
provides a lower limit [88, 93] on the axion decay constant, 

/ a >6xl0 8 GeV. (28) 

We identify the PQ scale with the axion decay constant. 6 Since our 
considered reheating temperatures are relatively large, it is probable that 
Tr. > Tpq ~ f a . 

If the PQ mechanism is supersymmetrised [94, 95], the axion a is part of 
a supermultiplet, the axion multiplet. It consists of the axino a containing 

6 This is equivalent to the choice N = 1, where N characterizes the colour anomaly of 
U(1)pq. With this choice we also avoid possible problems with topological defects. In the 
working scenario at the end the axion abundance from strings is negligible. For N 7^ 1 
our formulae hold, but f a = Jpq/N /pq. 
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the fermionic degrees of freedom, the axion itself and the saxino or saxion 
<f> sax , which is an additional real scalar degree of freedom. 

We investigate </> sax as candidate for the entropy-producing particle 4> of 
the previous sections. Thus the axion and the axino are unavoidable by- 
products, while it will become clear why they are no candidates themselves. 
However, the saxion is motivated from the strong CP problem and not from 
our scenario. 



5.2.1 Thermally Produced Multiplet 

After reheating, reactions like qq g4> sax and gg g^sax drive the sax- 
ion into thermal equilibrium if the reheating temperature is larger than its 
decoupling temperature [96] 

« = 10llGeV (ira) 2 O 3 - < 29 > 

where a s = gl{(j)/^K evaluated at the scale relevant for the processes under 
consideration. The equilibrium saxion abundance is given by 

^slx = 45C(3 j ^ 1.21 x 1(T 3 . (30) 
27rV(T s d a c x p ) 

Throughout this paper we assume for simplicity the particle content of the 
MSSM when we determine for example the relativistic degrees of freedom in 
the Universe. The numerical changes from adding the axion multiplet, for 
instance, are tiny as we can see in all estimates. 

The saxion becomes non-relativistic at a temperature T s ^ r x ~ 0.377n sax 

pfo 
nlsp 



around its mass [96]. From 0.37 m sax ~ T sax > T sax = X*? ~ m n i sp /25 



would arise a lower bound 

m sax > . (31) 

We will find that this is weaker than the lower bound on the saxion mass 
from early enough decay (36) and thus in nearly all cases and at least in the 
interesting ones does not yield any constraint. With (15) we see that if the 
saxion lives long enough, it dominates the energy density of the Universe 
below the temperature 

T= x = ~n e a> S ax * 1-6 GeV (^) . (32) 

We avoid matter domination during NLSP freeze-out by requiring T s ^ x < 
T^° sp ~ m n i sp /25, which gives an upper bound 

"••« <2 - 5TeV (l(j?5ev)- < 33 > 



19 



As we know (17), considerable entropy is produced only, if T™£ <C T sax . 
The saxion decay temperature can be derived from (13) with T = and 

^99 



t = l/rtf x , where 



2 3 

atmi 



rgg ^ s sax (<iA\ 

sax ~ 128tt 3 /2 1 j 

is the width of the dominant saxion decay into two gluons [96]. This 
yields [97] 

US g 53 MeV ( 10 ", GeV ) (lUm-) 1 (» ) f "ffA ' . (35) 



Here, a s has to be evaluated at m sax . As we do not consider an extremely 
large range of saxion masses, a s (m sa , x ) does not vary significantly. Besides, 
in the range of parameters considered, g*{T^£) remains approximately con- 
stant. Therefore, we drop the explicit dependence on a s and g*(T^) in the 
following equations. Together with the bound (21) from early enough decay, 
we obtain the lower limit 



T^after \ 3 / f 
min \ / J a 



2 2 

3 



^>™ &v liS Ural • (36 > 



If we compare this lower bound with (31), we see that (36) is stronger as 
long as 

2 

fa \ 3 > ( "tnlsp 

1.5xl0 10 GeV/ ~V10 2 GeV, 



(37) 



In any case, the saxion mass is constrained to a window. Since there is no 
additional source of SUSY breaking, one expects m sax ~ m susy and thus a 
saxion mass in the TeV range, i.e. 10 2 GeV < m SSLX < 1 TeV. Thus, from 
this discussion one might conclude that the requirements i)— iv) of Tab. 1 
are naturally fulfilled. 

Plugging (32) and (35) into (17) we obtain 



For simplicity we replace {gl^) for the moment by 2.2, the value estimated 
for (22). We plug in the bounds on the saxion mass (33) and (36) to find 

u ( fa \ ( W 2 GeV y ( f a y (AMeVy 

14 Vl0 12 GeVA ^msp J 0t> \W^GeVj \ J ' 

(39) 

The lower bound on A shows that (36) is always stronger than (31), since 
the inequality (37) is always true as long as significant entropy is produced 
and if this were not the case, (31) would not be considered at all. 
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If saxions are part of the particle spectrum, (39) shows two things, i) It 
is likely that saxions produce significant entropy in their decays. To avoid it, 
one would have to restrict the reheating temperature such that they never 
enter equilibrium, or to choose safe values for f a and m sax , e.g. m sax = 1 TeV 
and f a = 10 10 GeV. ii) The corresponding dilution factor is much smaller 
than the maximal value allowed by cosmology and preferred as a solution of 
the NLSP decay problem. 

The dilution factor can be increased by a larger axion decay constant, 
which makes the saxion more weakly interacting. From (39) we see that we 
need f a — 5.2 x 10 13 GeV to reach the maximum A ~ 0.75 x 10 3 of (22). 
This increases the decoupling temperature (29) and thereby the reheating 
temperature required to have the saxions in thermal equilibrium. If they 
did not enter equilibrium, the abundance would be y sax <C and the 

saxion would be useless for our purpose. From the requirement Tr > T s l x p 
and (29) we derive the upper bound 

f a < 1.0 x 10 12 GeV ( T ^-) " (—) 1 , (40) 

J a ~ V 4xl 12 GeVy V0.03/ v ; 

where a s (4 x 10 12 GeV) ~ 0.03. Already such a Tr corresponds — at least 
in the case of heavy gravitinos — to an allowed but relatively large dilution 
factor A ~ 10 3 , cf. (7) and (10). For the small As of (39) the situation 
becomes worse and is in fact inconsistent with itself. 

To sum up, if thermally produced saxions are to deliver the desired 
entropy, we need a large axion decay constant. Then we also need a large 
reheating temperature to make the saxion enter thermal equilibrium. This 
results in an overproduction of gravitinos, if they are produced without 
entering equilibrium, so that the scenario is not viable. 

On the other hand, if the gravitino is so light that it enters equilib- 
rium after reheating, the relic gravitino density becomes independent of 
the reheating temperature as mentioned in Sec. 3. Moreover, there could 
be another saxion production mechanism, which is the alternative we will 
concentrate on in the next section. 

Let us therefore continue discussing the requirements of Tab. 1, turning 
to the decay products of the saxion. Due to R-parity conservation, it must 
produce sparticles in pairs and thus cannot decay into single gravitinos. 
Besides, the decay into gravitino pairs is negligible, since it is suppressed 
by an additional factor of Mpj. Consequently, requirement vi) is satisfied 
without any effort. 

To fulfill requirement v) the decay into any other sparticle pair must be 
kinematically forbidden, i.e. m sax < 2m n i sp . This is the case if the saxion is 
lighter or not much heavier than the gravitino. Given that one expects both 
the gravitino and the saxion mass to be of order m susy , such a spectrum 
does not seem unlikely. It is understood that requirement v) does not apply 
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for a light gravitino produced in thermal equilibrium, since the NLSP decays 
early enough before BBN and does not overproduce gravitinos. 

On the other hand, the saxion may decay as well into two axions with [98] 

2 3 

paa ^ sax (»-• \ 

SaX ~ 327T/2 ' 1 j 

where the self-coupling x can be of order 1. In that case the Universe would 
be filled with relativistic axions during the process of entropy production. 
This would change the effective number of neutrino species, which could 
spoil the success of BBN. This requires i C 1 and there are concrete 
models realizing this [99]. 

Up to now, we went through the requirements i)-vi) of Tab. 1. We 
do not see any incompatibilities between the saxion producing entropy and 
gravitino dark matter. Hence, vii) is fulfilled automatically as well. Facing 
viii) we have to take care of the unavoidable by-products. 

Axion Due to the similar coupling strength the axion also enters equilib- 
rium, if the saxion does. Then its thermal abundance is the same as the 
saxion abundance (30). Due to its tiny mass [73] 

(1 10 CVV\ 
- j (42) 

its thermal relic density is negligibly small. The axion decays harmlessly 
into two photons with a lifetime many orders of magnitude larger than the 
age of the Universe [73] , while the gravitino cannot decay into axions due 
to R-parity. 

Axions are also produced by vacuum misalignment, which leads to the 
bound f a < 10 12 GeV to avoid overproduction [73]. However, for consider- 
able entropy production by the saxion this bound no longer holds, since the 
Universe is dominated by the saxion at the onset of axion oscillations. Then 
the axion density is given by [100] 

fi "^ a21 (ill)(ire) 2 - < 43 > 

If we require £l a /Qi)M = r < 1, we find 

h X<(JL\(^L\ (44) 
10 14 GeVy ~ V0.02/ \ J ' y J 

so that values of /„ > 10 12 GeV are indeed allowed. The bound from Q a <C 
^3/2 — ^dm is self-consistently cured by the decaying saxion. 

Altogether, there is no problem at all with the axion in our scenario. 
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Axino This is different for the axino a. Due to SUSY the axino has 
couplings similar to the saxion couplings and thus the same decoupling tem- 
perature (29) as the saxion and a similar relic abundance. The natural mass 
range for the axino is C?(keV) < < 0(m 3 / 2 ) ~ m, susy [101, 98]. We 
demand > m 3 / 2 to keep the gravitino as LSP. Then the mass range 
for the axino becomes similar to that of the saxion. In comparison to the 
saxion it has the opposite R-eigenvalue and thus must produce sparticles 
in its decays. With a light gravitino and axino NLSP, one would arrive 
at another NLSP decay problem. The situation would be worse than our 
starting point. If the axino should not produce gravitinos, which would lead 
to Y3/2 ~ y~ q = ^3/2) there must be another decay channel kinematically 
open, i.e. mj; > m n i sp . In the most interesting case m n i sp is close to the 
gravitino mass and the axino fulfills requirement vi). 

If kinematically allowed, the axino is expected to decay dominantly into 
a gluino-gluon pair with 

2 3 

From the requirement of early enough decay we derive the same lower bound 
on rria as found for the saxion (36). This would allow an axino lighter than 
the expected gluino mass. Weaker decays into sparticles and SM particles 
were to investigate. However, all these processes finally produce the lightest 
ordinary supersymmetric particle. The case of a heavy axino that decays 
after NLSP freeze-out has been studied in [102] for a neutralino dark matter 
scenario including the weaker decay into a neutralino and the re-annihilation 
of neutralinos. Even if the axinos are not in thermal equilibrium after infla- 
tion, they are regenerated by thermal scatterings and decays in the thermal 
plasma. The density produced in this way can be estimated in units of 
today's critical density as [103, 104, 105] 

The resulting neutralino density is many magnitudes larger than the thermal 
relic abundance (cf. e.g. Fig. 1), which in our scenario reintroduces the NLSP 
decay problem. In fact, the problem becomes much worse. Thus requirement 
v) is badly violated by the axino. 

If we require the axino to decay before NLSP freeze-out, so that we do 
not have to care about the produced number of NLSPs since they thermalize 
normally, we find by a derivation analogous to that of (36) the lower bound 
on the axino mass 

2 2 / , j , \ J 

6 



~ VlO^GeV/ Uo^GevJ WJ 



(47) 
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Since the gravitino problem could also be solved by making the gravitino 
comparably unnaturally heavy, such a large axino mass is not considered 
as a solution here. Furthermore, such an axino would produce considerable 
entropy with A^ ~ 290. From the discussion at the end of Sec. 3.2 we know 
that this would spoil our scenario, since A^ would dilute the saxion but not 
the NLSP. Thus, the situation is also inconsistent. The required axino mass 
to achieve A„ = 1 would be larger than about 10 7 TeV. 

Altogether, requirement viii) of Tab. 1 is badly violated by the axino. 
Consequently, the thermally produced saxion — and obviously also the axino 
itself — is ruled out as viable particle to produce significant entropy after 
NLSP freeze-out. The exception to this conclusion is a light gravitino in 
thermal equilibrium after reheating, since it allows for high reheating tem- 
peratures and the NLSP decay problem is absent. 

One may worry then if the strong CP problem can be solved by the 
Peccei-Quinn mechanism in scenarios of standard thermal leptogenesis with 
very light gravitino dark matter only. Going through the equations, espe- 
cially from (47) , we see that the axino becomes harmless for smaller axion de- 
cay constants f a < 10 10 GeV with an acceptable axino mass mj > 1.2 TeV. 
Since its decay into the gravitino is suppressed like (f a /Mp\) 2 , the contri- 
bution to the gravitino density from axino decay is negligible. However, by 
inspection of (38) we see that in this case the saxion is unable to produce 
a significant amount of entropy. Then also the axion abundance restricts f a 
to values smaller than about 10 10 GeV. 

In summary, by making the axino harmless we find that the thermally 
produced multiplet may also exist in scenarios of thermal leptogenesis with 
gravitino dark matter that does not enter equilibrium after reheating. How- 
ever, the axion decay constant is restricted to a small window. Moreover, 
the thermally produced multiplet is in fact useless for our purpose. This is 
due to two generic features of the considered scenario: i) Superpartners have 
similar couplings and masses, ii) The same coupling — or at least couplings 
of the same strength — were responsible for production and late decay of the 
entropy-producing particle. 



5.2.2 Generic Thermally Produced Particle 

The negative result for the saxion can be generalized to other late-decaying 
particles that are produced in thermal equlibrium by processes controlled 
by the same coupling as the decay. As the simplest estimate, let us assume 
that the particle <j> under consideration couples to SM particles via non- 
renormalisable interactions suppressed by an energy scale A and that the 
rate of reactions keeping <p in thermal equilibrium at high temperatures can 
be written as 

r p ; od = * S . m 
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where x is a model-dependent, dimensionless quantity containing couplings 
and kinematical factors, for example. The freeze-out from thermal equilib- 
rium occurs for H ~ r prod , which yields the decoupling temperature 

^..( ^K') Y A2 (49) 

* \ 90 J xMpi x x 10 18 GeV ' 1 ' 

For the decay we estimate 

rf c = yj t, (50) 

where y contains model-dependent factors. Generically, we expect x < y, 
where kinematical factors and the relation between number density and 
temperature tend to lead to a somewhat smaller x. For instance, for the 
saxion we find x ~ 6 x 1CP 7 and y ~ 3 x 10 -6 . We obtain the temperature 
after the decay as discussed in Sec. 3.2, 



/ 45 \ 4 (vm^M P ])2 „«2m? GeV 2 

T t C ~ 9 f H s 7 ' ~ 1.1 X 10 9 ^ , (51) 

1 27r 2 5 *(Tf c ) J A A 



assuming a sufficiently late decay to yield g*{T^ ec ) = 10.75. Together 
with the analogon of (32), which holds for any thermally produced scalar, 
and (17), we find the dilution factor 

A~l.lxl0~ 12 r, (52) 

(ym^ GeV) 2 

estimating as before (gl^ 3 ) — 9* (T^ ec ). 

Now we can use (7), (49), f^ 2 < Odm and Tr > T^ cp to obtain a lower 
limit on A and thus a constraint on the model parameters, 



A> — ( A V ( M ~9^) \ 2 A00 GeV \ 
~ x Vl0 14 GeVy V!0 3 GeVy \ m 3/2 )' 



i I i M I ( T< ' \ \ 

(53) 



3/2 / 

Furthermore, (17), (21) and (32) yield an upper limit on A, which can be 
combined with (53), resulting in 

A 1^1 AO 3 GeV\ 



Plugging this bound into (52) yields the maximal dilution factor that can 
be realized with a thermally produced generic scalar, 

\ ^ ( ( GeV\ / m 3/2 n! 
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Using further combinations of (21) and (51)— (53), we find that this maximal 
dilution is reached for 



Thus, we conclude that the generalized scenario allows for the production 
of some entropy, but we do not expect a dilution factor large enough to 
solve the NLSP decay problem. In order to avoid this conclusion, we have 
to consider a situation where the mechanisms for production and decay are 
different, so that the decoupling temperature and the decay temperature are 
no longer connected. 



5.2.3 $ sax as Oscillating Scalar 

As the saxion corresponds to a flat direction of the scalar potential lifted 
by SUSY breaking effects, it can develop a large field value during inflation. 
It begins to oscillate around the potential minimum when the Hubble pa- 
rameter becomes comparable to the saxion mass. This corresponds to the 
production of non-relativistic particles. The temperature at the onset of 
oscillations is 

r = .2.2*10.°Gev(^(^|)\ (58) 

Since we consider reheating temperatures higher than T s ° s ^ to enable thermal 
leptogenesis, the produced saxion abundance is independent of Tr and given 
by [106] 

osc i 
r sax _ rposc 



8 sax V Mpi 



where (ft sax denotes the initial amplitude of the oscillations and where we 
have assumed the simplest saxion potential, V = ^w-sax'/'sax- 

In this way production and decay are disconnected, circumventing the 
second feature mentioned at the end of the previous section. There is an 
additional free parameter, </>s ax - The saxion density is constrained by re- 
quirement iv), i.e. that it should not dominate before NLSP freeze-out. For 
the limiting case of domination onset at 7^° we obtain from (15) 

osc o o 

= ~ 4 T= X = Jzftp . (60) 
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Equalizing (59) and (60) we find for the initial amplitude 

( Of*) =6^ (61) 

\M Pl J r°- 1 ; 

or equivalently with T= x = T*j> sp ~ m n i sp /25 

(4Lc\ „,n^, m *( w10 GeV \ / 8.4GeV \*/ m nlsp xf ( g*(T^) \* 
V fa )~ 2 - txW { f a ){ ?Ti sax ) GeV J V 228.75 J ' 

(62) 

The easiest expectation for the initial amplitude is 0g ax ~ Mpi or <^ ax ~ f a . 
Interestingly, the estimate (62) yields an initial amplitude f a < 0g ax ~ 
\/f a Mp\ < Mpi, if we choose the harmless value f a = 10 10 GeV found above. 
According to (17) and (35), the maximal dilution (22) is achieved for a saxion 
mass m sax = 8.4 GeV on the lower boundary from early enough decay (36). 
The axion multiplet comes into thermal equilibrium after reheating, which 
gives the known limit > 1.2 TeV (47), avoiding problems with the axino. 
Thus, we have identified a working scenario where p° ax 3> psax, which enables 
significant entropy production while satisfying all requirements. 

Smaller f a are possible, too, provided that they respect the experimental 
bound (28). Larger f a and (/>g ax were not only in conflict with the scenario 
presented but also with standard cosmology. Furthermore, larger m sax are 
allowed, while they lead following (35) to smaller As. From the naturalness 
point of view, the required small saxion mass — compared to ttv^ and 7Tz SUS y — 
for maximal A might be the biggest concern. Nevertheless, we can conclude 
that the saxion as oscillating scalar can produce the desired entropy to soften 
the NLSP decay problem without violating any constraint from cosmology 
or observations. 

We would like to stress that our scenario does not contain more re- 
quirements than the standard scenario with axion multiplet but no entropy 
production. Instead, we only have to change the allowed windows for some 
parameters, most importantly ^J, ax and m sax . Avoiding axion overproduc- 
tion by vacuum misalignment even becomes easier. Other restrictions, in 
particular those on x in (41) as well as on /„ and mg, are the same as in 
the standard scenario. Note also that the initial amplitude of the saxion 
oscillations would have to be restricted to small values around f a , if one 
wanted to prevent entropy production. In other words, the classical field 
value of the saxion endangers the standard scenario. 

6 Conclusions 

We have discussed the possibility to solve the gravitino problem in scenarios 
with standard thermal leptogenesis and thus a high reheating temperature 
by late-time entropy production. Our setup has been stable gravitino dark 
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matter with a neutralino NLSP. We have estimated that thermal leptogen- 
esis is compatible with entropy production diluting the baryon asymmetry 
as well as the LSP and NLSP relic densities by up to three to four orders of 
magnitude. This amount of dilution roughly coincides with the maximum 
amount obtainable for radiation domination at NLSP freeze-out. 

For a gravitino LSP with a mass of 100 GeV, which allows for a reheat- 
ing temperature suitable for thermal leptogenesis, we have shown that a 
neutralino NLSP which is not much heavier can be diluted sufficiently to be 
compatible with BBN, i.e. its decays do not cause changes of the primordial 
light element abundances that are excluded by observations. However, this 
is only possible if the lightest neutralino contains a large wino or Higgsino 
component, whereas a bino-like neutralino remains excluded. 

We have discussed the general requirements on the particle producing the 
desired entropy and found that it is severely constrained. On the other hand, 
in some sense all these requirements either have to be fulfilled by long-lived 
particles anyway or are generically fulfilled. As a specific example, we have 
discussed the saxion from the axion multiplet. We have found that sufficient 
entropy production is not possible for a thermally produced saxion, where 
the same couplings are relevant for its production and its decay. This is 
due to two conflicting requirements: on the one hand sufficient production 
requires sufficiently strong couplings, while on the other hand sufficiently 
late decay requires weak couplings, where later decay corresponds to more 
entropy production. In the considered case, the allowed parameter ranges 
fail to overlap. Using simple estimates, we have generalized this negative 
conclusion to generic thermally produced particles. Furthermore, we have 
encountered severe problems with the saxion's superpartner, the axino. 

As an alternative, we have considered saxion production in coherent 
oscillations, which is independent of the saxion coupling. In this 
relatively light saxion with a mass around 10 GeV is indeed able to satisfy 
all requirements. Thus, if the Peccei-Quinn mechanism solves the strong CP 
problem, the potentially dangerous saxion decays can in fact turn out as a 
fortune, solving the a priori unrelated gravitino problem. 
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